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Abstract 

A passive orthonomic system of PDEs defines a submanifold in the cor- 
responding jet manifold, coordinated by so called parametric derivatives. 
We restrict the total differential operators and the prolongation of an evo- 
lutionary vector field v to this submanifold. We show that the vanishing 
of their commutators is equivalent to v being a generalized symmetry of 
the system. 

1 Concerning the status of this preprint 

After writing this preprint, I learnt (from an anonymous referee) that the Lie- 
bracket criterion is valid for all systems, not only orthonomic ones. Vinogradov- 
style derivation can be found for instance in [5j Chapter 4, §3]. The 'if part of 
the criterion follows from §3.3, while lemma 3.6 in §3.4 gives the 'only if part. 

Therefore, the present result is not the most general one. Still the reader 
may appreciate its Van der Kamp-style derivation. The intrinsic differential 
operators we will define, can be useful in practical situations. 

2 The standard symmetry condition 

In the majority of cases where exact solutions of differential equations can be 
found, the underlying property is a (continuous) symmetry of the equation [161 
HOj . And, in the theory of integrable equations, the recognition and classification 
methods based on the existence of symmetries have been particular successful 

[am nn a d. 

A symmetry-group transforms one solution of an equation to another solu- 
tion of the same equation. Although this idea goes back to Sophus Lie, we refer 
to [10] for a good introduction to the subject, numerous examples, applications 
and references. And we quote: 'The great power of Lie group theory lies in the 
crucial observation that one can replace the complicated, nonlinear conditions 
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for the invariance of the solution set of an equation under the group transfor- 
mations by an equivalent linear condition of infinitesimal invariance under the 
corresponding infinitesimal generators of the group action' |10j . In this paper 
we provide a characterization of symmetries that is different from the standard 
one, generalizing a similar characterization in the special setting of ordinary 
differential equations [THl eq. (3.35)], and evolution equations [TQ1 Prop. 5.19] 
to the setting of passive orthonomic systems. 

The natural framework in which symmetries of differential equations are 
studied is the so called jet-manifold M . Coordinates on M consist of p indepen- 
dent variables Xi , q dependent variables u a and the derivatives of the dependent 
variables, which are denoted using multi-index notation, e.g. 



u 2 



1,0,3 



drdt 3 



when x = (r, s,t) and u — (v,w). A typical point z € M is z = (xi,u a ,uj ( ). 
The ring of smooth functions on M will be denoted A. To indicate functional 
dependence of / 6 A we simply write f{z). Thus the system A(z) = 0, A 6 A n 
is a system of n PDEs. 

The action of a Lie group is defined on the space of dependent and indepen- 
dent variables, and then prolonged to an action on the jet manifold. Likewise 
the infinitesimal generator of the symmetry group is obtained by prolongation 
from an infinitesimal vector field on the base manifold. It turns out that any 
symmetry has an evolutionary representative [10\ Prop. 5.5]. In terms of the 
total differential operators 

a,K 

the prolongation pr^ of an evolutionary vector field Vq — ^2 a Q OL d/du OL is 

p^ = E^4' (2) 

A simple computation shows that these derivations on A commute among each 
other, we have [Di,Dj] = 0, i,j = 1, . . . ,p and 

[A,pr ] -0, i=l,...,p, QeA. (3) 

In fact, up to a linear combination of translational fields d/dxi, evolutionary 
vector fields are uniquely determined by property ([3]), cf. 10, Lemma 5.12]. 

The condition of infinitesimal invariance, the standard symmetry condition, 
is QB Theorem 2.31] 

pr Q A = mod A (4) 

in which case vq, or the tuple Q € A n itself, is called a (generalized) symmetry 
of the system A = 0. The tuple Q is a trivial symmetry if Q = mod A, 
which defines an equivalence relation on the space of symmetries. In section 
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[3] we show this is well defined for passive orthonomic systems. We restrict 
Q to be a function Q € B on the sub-manifold of the jet-manifold defined 
by our system of PDEs. Although this is a more than standard procedure, 
restricting the derivations to act on this sub-manifold is not standard at all, 
except possibly in the settings of ODEs and evolution equations. This is, at 
least from a philosophical point of view, not fully satisfying. 

In section [4j for any passive orthonomic system of partial differential equa- 
tions, we define intrinsic total differential operators Di and an intrinsic prolon- 
gation ptg, which are derivations on the sub-space B. Subsequently we show 
that the vanishing of the Lie brackets 

[D 4 ,pr Q ] =0, i = l,...,p, Q&B 

is equivalent to Q being a symmetry. 

3 Passive orthonomic systems 

Restricting to the sub-manifold is particularly easy when dealing with ortho- 
nomic systems, in which case this amounts to using the equations as substitution 
rules. However, in general the order of substituting and differentiating does 
matter, one encounters integrability conditions. For example, for the system 
u x = X, u y = Y to be formally integrable we need D y X = D X Y. In general, 
a finite number of integrability conditions suffices to make the system formally 
integrable, in which case the system is called passive. 

The idea of a passive orthonomic system is the main idea behind Riquier's 
existence theorems [17| and the corresponding algorithms for solving systems 
of PDEs due to Janet [2] . Riquier- Janet theory extended the works of Cauchy 
and Kovalevskaya, it takes a prominent place in computer algebra applied to 
PDE theory [12J, and it has lead to important developments in polynomial 
elimination theory [T]. The passive orthonomic system was the predecessor of 
what is now called an involutive system of PDEs. For our purpose, the concept 
of involutivity does not play a role. We adopt a similar philosophy as in [6], and 
stick to the setting of passive orthonomic systems. In that paper an efficient 
algorithm is given by which any orthonomic system can be made passive by 
construction of a sufficient set integrability conditions free of redundancies [6] . 

Let N = {0, 1,2,.. .} and N q = {1,2,..., q}. We denote the i-th component 
of J € N p by Ji and addition in N p is denoted by concatenation. A set of 
basis vectors for N p is given by {1,2, ... ,p}, where ij = 1 when i = j and 
ij = otherwise. Thus, with J,K € W, we have (JK)i = J, t + Ki, and in 
particular (Kj)i equals Ki or, when j = i, Ki + 1. Also, when L = JK we 
write J = L/K. Since total differential operators commute we can define a 
multi-differential operator Dk as 

D K = D^D^.--D^, (5) 

and we have Dku j = u K . 
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Choose n points (i a , J a ) G N g xN p , with nonzero J a , a — I, . . . , n. A deriva- 
tive wL is called principal if there exist L G N p such that (J, K) — (i a , J a L) for 
some a. The remaining ones are called parametric. The set of all (j, if) such 
that u 3 K is parametric is denoted S, and the subspace of A consisting of smooth 
functions of the parametric derivatives is denoted B. 

We also choose a ranking < on N g x N p , that is, a total order relation which 
is positive: 

VL,(j,K)<(j,KL), 
and, compatible with differentiation: 

(i,J)< (j,K)^(i,JL) < (i, KL), 

cf. Gang. 

We consider systems of n partial differential equations, with a — 1, ... ,n, 
u% =P a , P a G B. (6) 

The system ((6]) will be written shortly A = 0, where A" = vij a — P a . We make 
the following assumptions: 

i) the P a only depend on u° K with (j, if) < (i a , J a ), and 

m) J a K) = (i 13 , J^L) D K P a = D L PP. 

Such systems are called passive orthonomic systems. Their crucial property is 
that for any Q £ A, there is a unique Q G B such that Q = Q mod A. This Q 
can be obtained from the following reduction algorithm. 

Algorithm 1 Input: Expression Q £ A. Output: Expression Q G B. 

•k if no principal derivative appears in Q then return Q 

•k let u J K be the <-highest principal derivative appearing in Q, and let a,L 
be such that j, K = i a , J a L 

•k substitute u J K = Dj,P a in Q and call the result R 

-k return R 

The algorithm terminates because the highest principal derivative of R, if it 
exists, is <-smaller than u 3 K , due to assumption i). And, a different choice of a 
wouldn't change the result because of assumption ii). 

The following lemma states that differentiation is compatible with the above 
reduction A — f B, cf. Theorem 4.8]. 

Lemma 2 For any Q G A we have 

DkQ = D K Q 

Proof: Using a modified version of Algorithm [1] we can write Q = Q + -R(A), 

where R is some differential function of A such that R(0) = 0. Clearly DkR(A) 
vanishes. □ 
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4 The intrinsic symmetry condition 



Definition 3 We define an intrinsic multi-differential operator Dk ■ B — > B by 

® K P = D^P, PeB 
From this definition and Lemma [5] we obtain the following properties. 

Proposition 4 Intrinsic multi- differential operators are compatible with con- 
catenation, DkDl = ^KL- 

Proof: We have Dk^lP = D k dTp = DrDlP = DrlP = ®klP- □ 



Corollary 5 We have the analogue of equation {5p, Dk — £>i 1 ®2 2 ' ' '®p P ■ 

Corollary 6 Intrinsic total differential operators commute, [Sj, !0j] = 0. 

We would like to have a more intrinsic characterization of Di, that is, without 
reference to any principal derivative or total differential operator. For L e W 
we denote Sl — {(a, K) : (a, KL) e S}, which is a subset of 5. From equation 
fl} and Definition [3] it follows that 

3 (k,L)eSj L (i<*,J<*M/j)£S\Sj J°>M/j 

We note that when (k,L) G <S \ <Sj there exist a e N n , M e N p such that 
(k, Lj) = (i a , J a M), and, for any (3 G N n , N e W such that (k, Lj) = (i? ', J^iV) 
we have i a — i 13 and DmP" = DjqPP . Due to Corollary O equation ([7]) provides 
a recursive schema for intrinsic total differentiation. 

Proposition 7 The recursive schema for Di is well defined. 

Proof: We show the schema terminates using transfinite induction. For any 
Q e A, to evaluate DjQ, apart from some differentiation and multiplications, 
we need to evaluate a finite number of expressions DjP a . We assume that 
D^P 13 can be evaluated for all (i 13 , J&L) < (i a , J a I). Suppose P a depends on 
u 3 K . That implies (j, K) < (i a ,J a ). Suppose there are (3 G N n and L eW such 
that (j,KI) = (iP,jPL). Then, to evaluate £>/P a one may need to evaluate 
DlPP . By the induction hypothesis this can be done. □ 



Definition 8 We define intrinsic prolongation, denoted prg : B —> B, of an 

evolutionary vector field vq with Q G B by 

pv Q P = ^~P, PeB. 
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From equation @ and Definition [5] we get the intrinsic formula 

pr Q = J2 VkV^T- 

We now state and prove our main theorem. 
Theorem 9 A tuple Q G B is a symmetry of equation iff 

for all j . 
Proof: 

<= We calculate pigA" modulo A 

pr^A> = Dj^Qi- - pr^P" 

= D, /Q Q ia -pr Q P« (8) 

Next we calculate the commutator pDj,ptg]. Neglecting second order 
derivatives, we get 

%p*q= E ^o*^-' 

(fc,£)es L 



and 



(fe,i)es 3 - aUi (i«,j«M/i)e5\5^ ou J°<M/j 

Hence we get 

(i«,J°M/j)eS\S3 J<*M/j 

Suppose that J" 7^ 0. Then the action of the above vector field on u % Ja ,j 
yields the right hand side of equation ijHJ). Since we have chosen J" / £ 
N p this proves our case. 

Suppose Dj=Q iQ =px Q P a . Then 



(."J-M/jlES^- du J"M/j 



We will prove that [£ M ,pt Q ] P Q = for all a G N„ and M G N p . The 
statement is certainly true for A/ = 0. Assume that [Siv^ptg] P@ = 
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for all (iP,jPN) < (i a ,J a M). When M ^ there exists j such that 
M/j G N p . We write 

[2) M ,pr Q ] P a = J) j [S M/j ,pt Q ] P« + [%,pt g ] 2) M/j P Q 

The first term is zero by the induction hypothesis, so we concentrate on 
the second, which is 

Suppose ?" depends on u\. Then (fc,L) < (i a ,J a ). The function 
®M/jP a may depend on the derivative u k LM ^~ (namely, if u k hM ^ G <S, 
otherwise it depends on smaller derivatives). But when {i^,J^N/j) < 
(k,LM/j) < (i a ,J a M/j) by the induction hypothesis [J>N,pt Q ] P fi van- 
ishes. □ 



5 Discussion 

We have show that the infinitesimal condition [S3j,ptg] = 0, for all j, is 
completely equivalent to the more standard infinitesimal symmetry condition 
prgA = mod A in the setting of passive orthonomic systems ([6]). These 
conditions are sufficient for the corresponding group of transformations to be 
a symmetry, in the sense of transforming solutions into solutions [101 Theorem 
2.31]. In general, one needs to make non-degeneracy assumptions on the systems 
A to ensure that they are also necessary, see [TU1 Section 2.6 and 5.1]. However, 
in the setting of passive orthonomic systems local solvability is ensured by the 
existence theorems of Riquier and Janet |17[ [2] . This means that infinitesimal 
methods can be used for classifying all the symmetries of passive orthonomic 
systems. 

We do not propose that the bracket condition should replace the prolon- 
gation condition for practical calculations. On the other hand, the intrinsic 
differential operators do provide a convenient way of restricting expressions to 
the sub-manifold defined by the equation. For example, consider the Boussinesq 
equation 

l^tt — Pi P ^xxxx '^'U'x'U'XX 

If we let the prolongation prg, with Q = u txx — utu x , act on it we find 

Wq{uu-P) = Utttxx ~ tu t U 2 xx - 6u tx U x U xx + QutxxxUxx - Zutxxul 
^UtU x U xxx -\- Sllt xxxx U x -\- Qut xx 1l xxx -\- ^Ut x U xxxx 
~\~Ut^xxxxx ^txxxxxx ^^tx^tt IHtx'U't ^Htt^X' 

To see whether this vanishes modulo the equation we have to substitute iter- 
atively u u — P whenever we see two ^derivatives. An alternative is to work with 
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intrinsic differential operators. The parametric derivatives are u, u x , Ut Xl u xx , . . . , 
u x n,,Utx n i ■ ■ ■■ Clearly we have D x = D x . For the intrinsic total t-derivative we 
have 

s * = i^+y, 11 ^* — — 

at ou x n outx» 

n 

which can straightforwardly be implement in any symbolic manipulation pack- 
age. Without having to perform any further simplifications it follows that 

D?Q-pr Q P = 0. 

We do not want to check that pDt,ptg] = using a computer since it involves 
infinite vector-fields, which we can only deal with if we let them act on finite 
expressions. 

The infinitesimal symmetry Q is the first of an infinite hierarchy. The next 
member of this hierarchy is given by 

R "Z^xxxxx ~l~ 3u xx -\- ^U XXX U X TT^rEI 

5 o 

In fact, the symmetries of this hierarchy commute modulo the Boussinesq equa- 
tion, we have [prg,pr fl ] = 0. Using our bracket formulation of infinitesimal 
symmetry together with the Jacobi-identity one can conclude that the bracket 
of two symmetries should again be a symmetry. However, it is not a priori clear 
that this bracket is the intrinsic prolongation of something. 

We have observed that up to some scaling, the 'spacial part' of R equals 
the right hand side of the integrable potential Kaup-Kupershmidt equation [8 , 
Equation 4.2.5] 

15 2 20 3 

l^xxxxx l-Vll x U xxx "T ~^~^ xx "T ~^^xi 

and the spatial parts of the other symmetries in the hierarchy of the Boussinesq 
equation are the symmetries of the potential Kaup-Kupershmidt equation. At 
orders where the Kaup-Kupershmidt equation has gaps in its hierarchy, the 
symmetries of Boussinesq do not have a spacial part. 

Since the Frechet derivative of the Boussinesq equation is self-adjoint, its 
symmetries are also co-symmetries of the equation. Therefore they are in some 
sense dual to the Boussinesq equation. In the discrete setting this kind of 
duality was first introduced in [llj . Considering the dual equation Ut xx — Utu x , 
amongst its symmetries we find the hierarchy of the Sawada-Kotera equation 



U'l U'xxxxx dV'Xxx'U'X t ^ 



3 



3 

We note that a similar observation, using Lax-pairs, was made in [3]. These 
examples show the importance of extending the techniques developed for es- 
tablishing the integrability of evolution equations and classifying them, see 
[8j [14j [TBI HI El an d references in there, to the realm of passive orthonomic 
systems, or beyond. 
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